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INTRODUCTION 


This  report  contains  information  in  four  areas: 

a)  ^ 

^ /Air-Launched  Balloon  System  (ALBS):  Programs,  their  description, 
flowcharts,  and  Documentations* 

( y/ 

4 Balloon  Tracking  Using  VOR:  Description  and  Logic  of  M8themr>'cal 
Development,  and  a Program  Employing  These  Techniques  (Including 


^ ^ Flowchart  and  Documentation)^ 

^1~^  i Tethered  Balloon  Experiment;  Analysis  and  Recommendations^ 

-''llalloon  Plight  Simulations  Using  b(i/3in  AFB,  HOLLOMAN,  and 
WHITE  SANDS  Wind  Data. 
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I.  ALBS  Program  Documentation 


I n t roduct ion 

This  writo-up  describes  pi'ograms  pertaining  to  the  Air- 
Launched  Balloon  System  (ALBS),  which  were  originally  written 
for  a desk-top  calculator  (Wang  Model  452-1),  by  Andrew  Carten 
(AFGL) , ALBS  development  program  manager. 1 These  programs 
have  now  been  put  in  Fortran  so  that  they  may  be  run  on  the 
the  COC  6600,  and  this  section  provides  all  information 
necessary  to  run  these  programs.  For  each  program  are 
provided  a description,  listing,  flow  chart,  sample  deck 
set-up,  sample  output,  and  definitions  of  all  nomenclature, 
i.e.,  the  variabltj  names,  used  in  the  program.  As  much  as 
possible,  the  original  program  names,  subroutine  names,  and 
variable  names  have  been  retained.  The  underlying  mathematics 
and  logic,  unless  different,  are  not  described  in  any  depth 
here.  For  a more  comprehensive  treatment  the  reader  should 
consult:  the  Carten  report.^ 

The  first  program,  called  CHUTE,  treats  the  ALBS  system 
as  the  drogue  and  main  parachutes  are  being  deployed.  The 
second  program,  called  P14B,  treats  the  system  as  the 
balloon  is  being  inflated. 

Note  that  within  the  descriptions,  flow  charts,  etc. 
variable  names  are  often  referred  to  without  explanation. 

The  reader  is  referred  to  the  appropriate  section  on  'Nomen- 
clature ' . 
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Program  CMl/ni 


This  program  computes  thooret i ca 1 - performance  data 
(velocity,  drag,  dynamic  pressure,  etc.)  for  the  ALBS  system 
starting  from  the  point  when  either  the  drogue  parachute  or 
the  main  parachuco  is  beginning  to  open  until  it  is  fii^lly 
opened  (Oiti,TF)  . The  computations  are  based  on  an  initially 


chosen  opening  time  (TF"TFI).  The  volume  of  the  filled 
parachute  is  then  computed.  If  the  calculated  volume,  using 
the  initial  opening  time,  is  sufficiently  close  (less  than 
200  cu  ft)  to  the  theoretical  volume,  the  program  goes  on. 

If  not,  a now  opening  time  is  chosen,  and  the  first  calculation 
1s  repeated.  ^ 

When  an  acceptable  opening  time  has  been  found,  the  be- 
havior of  the  system  at  t^TF  is  computed.  The  calculations 
are  continued  until  the  velocity  of  the  system  is  'sufficiently 
close'  to  the  equilibrium  velocity  (VBHI).  The  criterion 
for  'sufficiently  close'  is  satisfied  when  the  change  in  the 
dynamic  pressure  (Q)  becomes  'very  small'.  The  system  data 
at  this  point  become  the  starting  data  for  another  program 
not  described  here. 

Tile  program  then  reverts  back  to  the  time  period  during 
which  the  parachute  was  in  the  process  of  opening  (O^t^TF) , 
ond  it  culculiitos  deceleration  and  opening  shock  forces. 
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Dutn  input  to  the  program  is  by  moans  of  a READ  statement. 
Each  time  the  READ  statement  is  executed,  l.e.,  each  line 
(or  card)  of  input  data  is  read  a new  set  of  calculations  is 
done.  Any  number  of  sots  of  calculations  may  be  done  at 
one  time.  In  this  way,  parametric  studies  of  any  of  the 
parameters  in  the  READ  statement  can  be  done  easily.  The 
data  must  be  'unformatted',  l.e.,  the  variables  must  be 
separated  by  commas. 

Each  time  a new  opening  time  (TP)  is  tried,  subroutine 
P4U  is  called  with  its  first  argument  equal  to  0,  and  one 
page  of  output  is  written.  When  an  appropriate  opening 
time  IS  found,  P4U  is  called  again,  but  this  time  with  its 
first  argument  equal  to  1,  and  calculations  are  made  for 
t>TF.  This  output  is  printed  on  the  same  page.  Subroutine 
PIO  is  then  called  to  calculate  the  acceleration  and  opening 
shock  forces,  and  this  output  is  written  on  a new  page. 

The  output  of  the  program  is  written  on  tape  1 and  tape  2. 
Tape  1 contains  the  output  of  main  interest,  while  on  tape  2 
is  written  supplementary  information  as  well  as  results  of 
calculations  used  for  checking  purposes, 

(Note  that  in  the  program  there  are  statements  with  ’CHECK' 
written  In  tho  first  5 columns.  Because  there  is  a 'C  in 
column  1,  they  are  treated  by  the  computer  as  comment  cards. 

In  order  to  use  those  statements  the  'CHECK''  must  be  deleted.) 
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PROGXAi  C-iUTE 


Tknk  0I>T«1 


^TN 


Qint/fT 


PRO;«*«  CHUTE  (UPElf  T*PEE|lNPjr,T*PESatNPUM 

«LJ$  PPPACHUTE  OECLOVHeNT  CALCULATIONS 
N'AirTEN  PON  A GARTEN  i APGL 

NOPLAT  VESPRINt,  ENNANUEL  COLLEOEi  i/TT 

JQHNON  PVEL/  V(21t|  OVLMl) 

CQNHON  /CONST/  PItCatOOiTPiCOiN 
OAtA  PtiG0/3.1Al99t  it.Zf 

READ  DATA 

READ  (lt*>  KNAX,  HOi  TFI  , il,  CO,  00«  CP»  TO 
IP  (EOP(9).NE.OI  STOP 

WRITE  riCAOINGS  FOR  TAPEZ 
WRITE  (fiZai  UATCiroOATt 
WRITE  lEtZe)  W,DO,COiCP 
WRITE  <2iZRI 

initialize  required  wariarles 

p . 1.  t TP  « TPI 
DTP  • .1 

50  » (PI*D0»»E»/lt. 

WLNAX  ■ <2,/(3.PPI»»  S)  »*D0*»  J 

ATART  CAlCULATIONS.ONC  /ALUE  OP  TP  FOR  EACH  K3UNT 
DO  lOD  FOUNT  ■ tiNPAX 
call  PLJ  (0«HO)H|V0iWPt|d) 

VTHI  ■ SQRT<W/(,001ia9*SIGHA(HO)»CO*SGM 

AUNROUIlNES  P9  AND  PPU  IN  A.  GARTEN'S  PROGRAMS 
/L  • Q< 

:i  « oo*p»*TP/pi 
00  me  j I iiu 
r ■ FLOAT! 

9WL<J»  » -Cl*W(J>»m.-T> ‘T**  IR./S.I  - Z .•CP'T  • U , -T» » 

VL  ■ VI  A .1»DVL(J» 

SONTINUI 
WRITE  (ZiZei  V 
WRITE  (tiZ«l  DVL 

IS  TP  ACCEPTAOLET  IP  IT  tSi  GO  TO  ZOO. 

OIPPVL  • ARSIVLHAX  • VL» 

IP  (OIPPVL.LE.ZQOit  GO  TO  100 

P w IVLHAX  - VD/3IPFVL 

WRITE  (ItlCI  KOUHT,TP|VLHAX,VLt3tEPVL 

WRITE  (Z,10l  FOUNT, TP, VLNAXi  VLtOIPPVL 

IP  (DIFPVL.LE.IOOD.I  DTP  • • OZ 

TP  « IP  ♦ P»DTP 

:qnt  inu! 
lONTlNUE 

WAKE  CALCULATICNS  FOR  TIWE,GT,TP 
call  PLU  (1,H,HP,VF1,VP,AI 

WRITE  a,U)  TP,  VEHI,  VLHAX,  VL 
WRITE  (Z,I0)  FOUNT, TP, VLHAXiVL,3IPFVL 

WAKE  CALC  OP  ACC  AN(J  FORCES,  WHILE  CHUTE  IS  dEING  OEHLOYEO. 
CALL  PIO  (HOT 


PORWAT 
FORMAT 
l.//»ll 
PONMAT 
FORMAT 
l,»  STS 
Z,«  OOE 
FORMAT 
FORMAT 
STOP 
END 


(IHO,  t9,SPlO,Z> 

(1HD,///|UX,  *0PEN1N6  TIME  ■ * ,P  t,  Z , OX  ,•  V£H  X ■ *,Pr.Z 
X,wVOLHAX,  * *,Pia.Z,lCXi'CALC  VOLUME  ■ >,P10.Zt 
(IH1,*ALBS  ••  SUPPLENENTART  INFORMATION*,  I 0X«A  IQt 
(IHO,  10X,*tNPur  VAR1A9LIS  NOT  LISTED  ELSEWHERE*,// 

TEH  HEIGHT  • *|PS.1,1SX,'  OlAN  OP  CHUTE  • *,PB>3,// 

PP  OP  DRAG  * *,Pt.3,lJK,'  C OF  POROSITV  • •»FS,3,//) 
tlHO,  10X,*OTHER  GALCULAIED  VANlABLES*t 
( IHO,  llFl  0,Z> 


000100 

oootio 

000  IZD 
003  133 

ooomo 

000190 
000  ICO 
00D17D 
000100 
000190 

ooozoo 

OOOZIJ 
003ZZD 
OO0Z3O 
OOOZwO 
OOOZSO 
ODuzeo 
000  Z/0 
000100 
000190 
000300 
000310 
0003ZO 
000330 
ddC3l|3 
003  3 03 
3011313 
003373 
300300 
003 190 
OOOLOO 
000410 
J3C4Z3 
000430 
000440 
000400 
30C4H 
0004TO 
000400 
000490 
30CS33 
OODIIO 
Q009ZO 
003  A30 
OOOStO 
OOOIVO 
OOOOsO 
OOOOPO 
000903 
aoo>90 
oooooo 
000010 
ODOOZO 
00  3 0 33 
OODOWQ 
000610 
000460 
003070 
OOOOdO 
003010 
000700 
000710 
0007ZO 
300  730 
000  7 40 
000  790 
000/  00 
000770 
OOO7R0 
003790 
OOOAOO 
Dooaio 


IMlkhMAUtU^IOBtikr^  M^IIa  9 ^li 


SUCAOUTInE  Pio 


(4/74 


OPT*l 


rr-i  4,5»(«i4 


OS/U/77  1MZ.3S 


1 


! 


10 


20 


2S 


30 


31 


SU'i^OUriNE  PIO  <Ht  0U1S20 

COMMON  /tfEL/  V(21I|  OVLIllI  001S30 

common  /CONST/  PI,Oa,OOiIF,CO,(<  001S40 

OtMlNSION  CVOTUll  OOlkSO 

flOllOO 

CC  INITIALIZE  ftCOUlKEiJ  <t41tOLES  001170 

OOSOH  ■ CO*(PI/i*.»»00*»2  001580 

1 ■ (H*100000il/(Siart«(NI*CO*00**3>  001590 

9 • 1204»CDSON/00»*  3 001500 

: ■ 1*00  001510 

001510 

CC  M9ITE  HCAOINCS  ON  TlPKl  001530 

KIITE  (lilt  OATEUCOITI  001540 

001550 

SC  STMT  SttCULlTlONSi  INTEKVUS  OF  TF/10.  001640 

UO  100  I ■ lill  001670 

r • FLOlTd-D/lO.  001650 

DVOnil  • <22.9*V(II  - 0*TF*I*vm**2  ♦ C*TF»/(1  * 22,i*TI  001590 

F • -(M/ (iO*TF|)  •ovor  (1>  001700 

F0«  F » M 001713 

1 • FO/H  001720 

CC  M4ITE  CNE  LINE  OF  OUTPUT  00173J 

MRITE  UilO)  T,  OVKIlt  OVQTmi  F.  FO  • 3 001740 

too  COnTINUI  001750 

001750 

:e  M91TE  ICOITtONAl.  INFOAHiriCN  001770 

MmTE  u,20»  It  B,  C)  COSCM  001700 

001790 

1 F0.(HAT  (IHl,  10X,*P905HAm  P9*  « 1 3X,  •PkOCRAN  PIO*  . 2 0 X > AlOi  / 7//  00  10  00 

l,llt*.*T/TF  OVOt  tcu»n  07/0T  IFPS/ S»  *.  lOX.  *»  ( Ld)  * 1 9X  » 001010 

2 *FCILei*i7Xi*F0/M  001020 

10  FORMAT  (lH0il0XfFS(l»5FlU.2|F14,31  001030 

20  FORMAT  (1M0|10X|*A  • • • F II.  3 » 5X  i*8  it  »>F14.l|Mi*0  • *|F  12.  3 » 7/001 140 

1|3X.*CCS0  <•  *|F12.2I  0010  50 

4CTUFN  OOlOeO 

END  001570 


function  SliHl 


74/74  OPT«l 


FTN  4.5F..14  01/12/77 


1 


5 


10 


funCTILM  SISHA  (H)  OOlOOO 

DIMENSION  510(  251  001590 

THIS  FUNCT104  INTERPOLATES  lETHEEN  VALUES  OF  THE  DENSITY  RATIO,  001900 
OIVEN  NY  SIQ  IN  intervals  CF  ICOG  FT.  001910 

DATA  SID/. 97107,. 94275,. 9151  3, .A5011,. 55170,. 831)0  001)20 

1..  51  07  0,.  7 54  09,  .75205  m 73  059,  *71550,.  59333..  67 151  0019  30 

2..  05022..  51945..  51921,.  55  945,.  07021,. 591 441.93315  001)40 

3mH934,.49796.'.45105m45452.  .44519/  CODIj 

J > H/1000.  001940 

XJJ  ■ H 4 1000.»fluAT(  Jl  0019  70 

SIDMA  • SIC(J>  * ISlSUfll  • SIC(JI)*XJJ71D00.  001)50 

RETURN  00199} 

•NO  002000 


10 


13. 12.36 


...  i.  Ol'iCl.I.L-!. 


0PT»1 


FTN 


U«t2«34 


suaRouTr^g  p<»j 


fkffk 


1 

SUOKOUTINE  P4U  (IFL*S«HQiHF,VOiV'»SI 

001029 

cc 

IFLAG  I 0 — CHUTt  IS  3E1NG  OEPtOTED 

0009  30 

cc 

IFLtC  < 1 ••  CHUTE  HAS  dEEH  DEPLOYED 

000949 

COHHON  /CONST/  PI,Gll,O0<TF,CO»H 

900SS9 

s 

SOHPON  /WEL/  V(2U,  CVLdli 

909969 

000979 

INITIALIZE  REQUIRED  YARIA8LES 

900999 

30  ■ (FI*D0»*2»/<.. 

0009  90 

1 « .O01189«SICH*(H0l*V0**2 

OOOvOO 

10 

oobno 

WRITE  hEtOINSS  FOR  TAPElt  Ir  NECESSARY. 

000920 

IF  IIFLAG.EQ.il  60  TO  SO 

0009  30 

MRITF  (1,101  OATEIT03AYV 

009|il0 

WRITE  (1,21  R,  Q,  YO,  HO 

000  9 93 

13 

30 

CONTINUE 

000990 

000970 

Y(  10*IFLA6a1I  • VI  • VI  ■ VO 

000390 

N I HO 

000990 

31  > C. 

001030 

20 

UT  « TF/50. 

QOIOID 

OVl  ■ CVO  ■ *G0«OT»(l.  - Q*C3*SO»FLOAT  (IFLAGI/MI 

901020 

I9d  ■ SO'IFLAG  « 1 i IFF  • IB3  * 49 

D01Q33 

0010  40 

CC 

START  calculations,  INTERVALS  OF  OT  ■ TF/50. 

001090 

25 

DO  100  I ■ ISBtIFF 

001060 

r > FLOATdl/SO. 

001  3 73 

TINE  ■ T»TF 

001090 

A ■ (SO'OI  •AHINKT.I.I  * a 

001090 

3VPRE  1 2. ‘OVl  - OVO 

OOllDO 

30 

VPRE  • VI  ♦ OVPRE 

OOlllO 

VEAV  * (VI  * VPREI/2. 

001120 

OHORE  ■ VEAV’OT 

001130 

3 ■ .00llfl9»SIGHA(W»0HPRE)»VPR£»»2 

001140 

•j  1 O*C0*A 

001190 

33 

:vi  ■ -Gfl»0T*(l.  - OFWI 

001190 

VI  ■ VI  A OVI 

001170 

3H1  ■ .!*(Vt  ♦ Vll’or 

OOllAO 

W • H ♦ OHI 

OOI19V 

VI  «VI 

001209 

LQ 

3V0  > CVl 

001210 

jvi  » CVl 

001220 

1)01233 

:c 

WRITE  ONE  LINE  OF  OUTPUT,  EVERY  FIFTH  CALCULATION 

001249 

II  < MC0(1,9I 

001290 

kS 

IF  (It.NE.  01  GO  TO  90 

001290 

IJ  ■ t/S  A 1 

0012  70 

VdJI  ■ VI 

C0L20C 

WRITE  (1, 1)A,T, TINE, Q, OVl, VI, U,H 

0012)0 

001300 

so 

cc 

FOR  IFlAG  ■ li  calc  STOPS  WHEN  OlFF  IN  3 IS  SNALL 

001310 

IF  (IFLAG. EO. 01  G3  TO  90 

031323 

OIFFQ  ■ AeS(Q  - QII 

001330 

IF  (01FFi].LE..a01l  RETURN 

001340 

QI  • Q 

001S50 

ri 

90 

CONTINUE 

001390 

0013  70 

-;h£C 

< IF  (IFLAG. EO. 11  WRITE  (2,11  A,  T,T1H£  ,3  ,0V  I , V 1 , 0 ,H 

001360 

CH£C 

K WRITE  (2,11  A,  T,  TINE,  Q,  OVI,  VI,  0,  H 

001390 

100 

CONTINUE 

001400 

f<0 

OOlxLO 

cc 

RF  ANL  VF  are  output  OF  P4L  SENT  SACK  TO  AL9S 

00l4i!0 

NF  > H ( VF  ■ VI 

0014  30 

WRITE  (2,21  S)0,VF,RF 

001440 

1 

FORMAT  (1H0,10X,F11.  2,F9.  t.F9,2,2Fll,4,3F11.2t 

001413 

r>y 

» 

FORMAT  (IH  ,lQX,Fll.2,l/K,F11.4,llK,2(F11.2,llVn 

0014LO 

to 

FORMAT  ( 1H1,10X,*PROGRAN  P4U  *,AtO,//// 

001470 

1,12X,»AREA(SQFT)  T/TF  T(SECI  Q(PSFI  OV(FPS) 

V (F*S)001460 

2 C(LB>  H(FT»*,//I 

301490 

RETURN 

00 1100 

20 

ENO 

001910 
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] 

I 
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-i 
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SUBRQUTtHE  R4U 


1 PiITIAL  CAlCULATlCNal 


■^HEADINGS 
POR  NEW 
«LCULftT10NS 


I WRITE  HEADINGS  FOR  TA^ 


MAKE  CALCULATIONS 
OEFIME  LIMITS  OF 
00  LOOP  TO  FOLLDVI 


mis  SUBROUTINE  NWKES  CALCULATIO^lS 
AS  ms  CmUTE  is  DESCENDING, 

BOm  WHILE  IT  IS  OPENING  AND 
AFTER  IT  ^^^S  OPENED. 


MAKE  CALCULATIONS 
.AT  INTERVALS  OF  . 
At  « TF/S&>^ 


a-IE  PASS  FOR  EACH  TIME. 


/•^SHOULD  ^ 
OUTPUT  BE 
s^PRlNTEOT^ 


WRITE  ONE  LIN!  OF  OUTPUT  EVERY 
FIFm  CALCULATION. 


PUT  CURRENT  VELOCITY 


WRITE  ONE  LINE  OF 
(output  ON  TAPE  1 


IS  PARACHUTE  FULLY  OPENED? 


CALCULATE 
CHPNGE  IN 


IS  VELOCITY  SUFFICIENTLY  CLOSE 
TO  EQUILIBRIUM  VELOCITY  (VtHI)? 


change  IN 
V SMALL>' 


SET  FINAL  HEIGHT  AND  VELOCITY. 


RETURN 


SUBROUTINE  PIQ 


THIS  SUBROl/riNE  COMPUTES  OPENING 
SHOCK  FORCES. 


ONE  PASS  FOR  EACH  TIME  INTERVAL. 
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■ end  of  file 
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Fr'.;  .-.HUTE;  aHPEIMAF'EE)  IMFUT^TAPES'IMFUT/ 


end  of  record 
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11111231 
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1121 

miimmii 

nil 

1 1 1 1 1 1 1 

IGlinil 

jniiim 

1 1 0 

00  4 

111  11  II 

2 n II  M n I 

iiniii 

1111 

1111111111111 
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2 2 2 2 2 2 J 

21222222 

2 2 1 2 2 2 2 2 2 
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mil  n 

1 1 m 3 m 1 

1 1 2 1 III 

nil 

1 1 1 1 1 1 1 1 1 1 1 1 1 

nil 

iiniii 

1 1 1 1 1 12  1 

IIMIIIII 

1 1 0 

1 1 1 

limn! 

1 1 1 1 1 1 1 1 

9232:12322 

t t II  1 1 <1  11  1 II  .) 
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n II  It  II 
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Nomenclature 


Projjrum  CHUTE 
B 

CD 

CP 

Cl 

DIFFVL 

DTF 

DVL 

DO 

F 


GO 

H 

HF 

HO 

J 

KM  AX 


Cross-sectional  area  (b)  of  parachute 
at  t«0 

Coefficient  of  drag,  Cq 
Coefficient  of  porosity,  Cp 
(Cl)  - d^^j/iT 

Absolute  difference  between  the  calculated 
and  theoretical  parachute  volume 

Correction  to  tf  to  get  DIFFVL  below  an 
arbitrary  limit 

Rate  of  change  (dV/dT)  of  volume  with 
T ■ t/t^ 

Parachute  diameter,  d^ 

0 

Number  used  to  decide  whether  t£  should  be 
Increased  or  decreased  to  produce  an 
acceptable  parachute  volume. 

Acceleration  of  gravity  (g^  ■ 32,2ft/scc^) 

Final  height  returned  from  subroutine  P4U 

Final  height  after  chute  has  achieved 
equilibrium  velocity 

Initial  height, 

Index  used  to  increment  successive 
iteration 

Maximum  number  of  different  opening  times 
(t^)  that  will  be  tried  for  a given  run. 
This  is  an  additional  safeguard  against 
the  program  going  into  infinite  loop. 


K^UNT 

PI 

SIC 

SO 

T 

TF 

TFI 

T0DAY 

V 

VEHI 

VF 

VFI 

VL 

VLMAX 

VO 

W 


Index  used  to  count  the  number  of  different 
opening  times  (t^)  that  have  been  tried. 

TT 

Array  defining  a at  1000-ft  intervals,  whore 
SIG(I)  « O(h-1000*lj,  up  to  25000  ft. 

Cross>sectional  area  (S^)  of  parachute  when 
full/  opened  (t>tf) 

Fraction  of  opening  time,  T ■ t/t^ 

Opening  time,  tf 

Initial  value  of  opening  time,  read  into 
the  program 

Dummy  variable  used  to  print  date 

Array  containing  values  of  velocity  (v) , 
where  V(l)  ■ v (T» , 1) , v ( 2) -v(T« . 2) , etc. 

Equilibrium  velocity  of  parachute 

system 

Final  velocity,  VFeVBHI 

Velocity  after  parachute  is  fully  opened 

(t"tf) 

Parachute  volume  calculated  using  opening 
time  t*tf 

Theoretical  parachute  volume 

Initial  velocity  of  parachute  system  (at 
t-0) 

Weight  of  parachute  system,  Ws 


I 
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I 4 

( 
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21 


Nomenc 1 aturo 


Subroutine  P4U 
A 


B 

CD 

D 

DHI 

DHPRE 

DIPFQ 


DT 

DVI 

DVL 

DVPRB 

DVO,  DVI 

DO,  GO 
H 

HF 

HO 

I 


Cross-soctional  area  of  parachute  at  a 
given  time  t 

Same  as  in  CHUTE 

Same  as  in  CHUTE 

Drag  of  parachute,  D 

Change  in  height  in  time  DT*tf/50 

Initial  calculation  of  change  in  height 

Change  in  q in  time  ta.ltf.  This  is  used 
onlx  after  the  chute  is  fully  opened 
(IFLAG  ■ 1)  to  determine  if  the  velocity 
CV)  is  close  enough  to  the  equilibrium 
velocity  (VEHI) 

Time  Increment/lteratloni  DT  ■ tf/SO 
Change  in  velocity  (V)  in  time  DT 
Same  as  in  CHUTE 

Initial  calculation  of  velocity  change 
in  time  DT 

Used  in  determining  DVPRE.  They  are  the 
two  previous  DVI's 

Same  as  in  CHUTE 

Height 

Final  height  when  P4U  returns  control  to 
CHUTE 

Initial  height  when  P4U  starts 

Index  used  to  Incremont  successive  iterations 
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IBB,  IFF 
I FLAG 


II 


IJ 


PI 


q 

Ql 

so,  T,  TF 
TIME 
T0DAY  ,V 
VEAV 
VF 


VI 

VPRE 


VU 

VI 

W 


Beginning  and  final  values  of  index  I 

Index  used  to  tell  P4U  whether  the  parachute 
is  in  the  process  of  openinii  (IFLAG  ■ Qj 
0_<t£tf5  or  is  fully  opened  (IFLAG  ■ 1,  t^t^) 

Index  to  determine  whether  a particular 
Iteration  should  be  printed.  (Only  1 out 
of  S is  printed) 

Index  used  to  place  a velocity  value  (VI) 
into  the  array  V;  this  index  corresponds 
to  a time  such  that  V(IJ)»v(T),  1J«10*T+1 

Same  as  in  CHUTE 

2 

Dynamic  pressure,  q ••  *j^v 

Previous  value  of  Q;  used  to  calculatu  DIFPQ 
Same  as  in  CHUTE 
Real  time  t(sec) 

Same  as  in  CHUTE 

Average  velocity  over  the  interval  DT 

Pinal  velocity  when  P4U  returns  control 
to  CHUTE 

New  velocity  at  end  of  time  increment  DT 

Preliminary  value  of  now  velocity  at  end 
of  time  increment  OT 

Initial  velocity  when  P4U  starts 

Used  to  store  the  previous  velocity  (VI) 

Same  as  in  CHUTE 
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r Trjsufjiirtytinftr? 


Nomonc 1 aturo 

Subroutine  PIO 
A , a , C 
CD 

CDSOM 
DVDT 
DVL,  DO 
F 

PO 

G 

H 

I 

PI  ,T,TIM|3DAY.V,W 


Coofflclcnts  used  to  calculate  DVDT 

Same  as  in  CHUTE 

(C„S  ) - C„S 

'•  D o-'m  D o 

Acceleration  of  parachute  system  (ft/sec^) 
Same  as  in  CHUTE 

Force  produced  by  the  drag  of  the 
parachute  (IbsJ 

Sum  of  the  drag  of  the  parachute  and  the 
force  of  gravity 

Acceleration  in  g's  corresponding  to  FO 

Height  at  which  a is  to  be  calculated 
This  value  o(H)  Is  held  constant 
throughout  the  subroutine 

Index  used  to  specify  the  time  (T".1*I) 

Same  as  in  CHUTE 


Function 

H 

J 

SIG 

SIGIIA 

XJJ 


SIGMA 


Height  at  which  SIGMA  is  to  be 
calculated 

Index  used  to  determine  1000-ft 
interval 

Same  as  in  CHUTE 

Value  of  o(H)  calculated  and  returned 
to  calling  program  or  subroutine 

Variable  which  is  used  to  interpolate 
0 between  h»1000J  and  h"1000(J+l) 


1 


A 

1 


1 

■/ 

I! 
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Program  P14B 


This  program  treats  the  ALBS  system  starting  from 
the  point  when  the  balloon  has  been  completely  extracted, 
and  the  helium  has  begun  to  enter  the  balloon.  Calculations 
are  done  every  200  ft,  as  the  system  descends  through  the 
atmosphere,  and  one  line  of  output  Is  printed  for  each 
height.  The  effect  of  the  system  weight,  the  drags  of 
the  drogue  and  main  parachutes,  the  lift  produced  by 
the  Increasing  volume  of  helium  In  the  balloon,  and  the 
change  in  temperature  and  pressure  of  the  atmosphere  are 
all  taken  into  account  by  the  program. 

Data  is  introduced  to  the  program  thru  a READ 
statement.  Bach  time  a card  of  data  is  read,  another  set 
of  calculations  is  done.  Any  number  of  sets  of  calculations 
may  be  done  during  the  execution  of  the  program.  In  this 
way,  parametric  studies  of  any  of  the  parameters  in  the 
READ  statement  can  be  done  easily.  The  data  must  be 
'unformatted',  1.  e.,  the  variables  must  be  separated  by 
commas . 

The  output  of  the  program  is  written  on  tapes  1 and  2 
(See  Figures  1-4).  After  the  headings  are  written,  one 
lino  of  output  is  written  for  each  height  C200-ft  intervals). 
Since  the  values  of  so  many  variables  are  written,  the 
width  of  one  page  of  computer  paper  is  not  enough  to  write 
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them  all  across  the  page  In  one  line.  Hence,  for  each 
height,  some  of  the  variables  are  written  on  tape  1 
(Figures  1 and  3)  and  the  rest  on  tape  2 (Figures  2 and  4). 
After  completion  of  the  job,  the  pages  are  separated  and 
each  set  of  calculations  is  taped  together  so  that  the 
data  for  a given  height  can  be  read  across  the  page  in 
one  line. 

Each  time  another  READ  statement  is  executed,  another 
set  of  calculations  is  done.  For  each  new  set  of  calcula- 
tions, the  headings  are  reprinted  at  the  tops  of  new  pages 
on  Capes  1 and  2. 

The  Fortran  version,  explained  here,  also  contains 
two  new  features  not  found  in  the  original  version  written 
for  the  desk-top  calculator.  The  first  takes  into  account 
the  adiabatic  heating  of  the  gas  in  the  balloon  as  the 
system  descends  into  a denser  and  denser  atmosphere.  The 
second  provides  an  interpolation  to  determine  the  exact 
iieight  at  which  the  helium  is  exhausted  (instead  of  the 
height  of  the  nearest  200-ft  interval)  and  the  values  of  the 
other  parameters  at  that  height.  These  features  are  explained 
in  detail  below. 
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Pp»05^4H  P14J 


?4/f4  OPT»l 


PTN  4.S*414 


09/12/77  13.12.49 


I 


9 


10 


19 


20 


29 


30 


39 


<■0 


41 


! 0 


95 


P90CRAH  P14B  (TAPCl t TAPE2 i INPUTiT APe9*lNPUT ) 

000100 

cc 

OOOlltl 

:c 

AL9I  balloon  INPLATION  CALCULATIONS 

000120 

oc 

HRITTCN  rOK  A.  CAftTCNt  AF6L 

000130 

cc 

ST  ROBERT  VCSPRXNIi  CHHANUEL  COLLECEi  2/77 

000140 

cc 

dODlIl 

CI9ENSI0N  PRCSKXfli  S13<2B» 

ODOlOO 

CATA  P4ESS/.9B43M92BA|tlBA2> .BASA* < B320 ». 0 014 • • 7 717 t . 7420i 

0001  70 

1 .7t49«.BI7B».BBlAiiB3eiM6119».BB7AMfiA46MB422i.l206i 

000160 

2 .<»997i>47«9i  .4B99i.4410i  .4227 M'ltBii. 30»0«. 3716/ 

00019) 

CAT  A SIS/.  97107,  .94270,  .91B1 3 , .tOBll, . OB  17  0 , . 1 399  0 

000200 

1,. 4 1070,. 76609,. 76200,  .73099,. 71900,. 693 93,.  67191 

000210 

2, . 69  022, . 62940, . 10921, . 90940,  .97021,. 99144, . 93310 

000220 

3,.  9 1 93  4,.  49  7 90,.  40  1 00,.  4 6 4 02,. 44  0 9 9/ 

J03230 

3ATA  PI, Cl, OH  /3. 14119, 3. 0413, -200./ 
4EA0  OATA 

000240 

• >4*4  »4 
4*W«0  W 

0002SO 

!0 

OEAO  IC,*t  HI,  MSRE7,  THE,  TR,  TREF,  90,  DT,  TLIEL,  HO,  DAD, 

OAH  000200 

IF  (E0F<9) .NE.OI  STOP 

COCO  )0 
000260 

oGs; 

INITIALI2E  REQUXREO  VARIaOLES 

000290 

T44C  • 102.44 

000300 

OIEHP  ■•01  TEH*  ■ THE 

C06S1C 

IFLA3  ■ 0 S F ■ 1. 

000320 

ITOP  ■ IFIX(HX/200.>  • 1 

000330 

IF  (ITOP. 6T. 124)  ITOP  « 124 

00(  340 

I • ITOP  ♦ 1 

000390 

J ■ X/l 

000360 

JJ  ■ I - 9*J 

Q0D370 

SICHA  ■ 81C(J)  * <SI3(Jtll-SI0(J))«FL0AT(JJ)/l. 

40C36Q 

P * PP.ESS(JI  * <PRESS<JT1)-PRESS(J))»FL0AT(JJ)/9. 

000390 

/CL  • TM  ■ T M 0. 

000400 

3AS  • OAO  * OAH 

000410 

000420 

MRITE  HEAOINCS  And  OATA  AT  INITIAL  HEX5HT 

000430 

WRITE  (1,10)  OATE(TOOAY) 

000440 

WRITE  (1,11)  THE, TREF, TRiOT 

000490 

WRITE  (2,12)  OAH,  OAO,  MU,  TLREL 

000460 

WRITE  (1,1) 

000470 

WRITE  (2,2) 

000400 

WRITE  (1,13)  HI,  SIGHA,  P,  TO 

000490 

WRITE  (2,141  HI,  WSREF,  OAS 

000900 

000910 

cc:: 

start  calculations,  20o<ft  intervals 

000620 

30  1000  11  ■ 1,IT0P 

000930 

I • ITOP  - 11  ♦ 1 

000340 

J » 1/9 

001990 

OH  » OT*TR  \ 

OOOIbO 

THM  • TH  ♦ DM 

000970 

000960 

cc:: 

HAS  HELIUM  BEEN  EXHAUST'OT  IF  NOT,  50  T 0 60. 

000990 

IF  (TMM.LT.TmHE)  30  TO  80 

000600 

IFLAC  ■ 1 

ODOblO 

F ■ (TMHE«TM)/0M 

00062) 

T.H  • THHE 

000639 

SO  To  90 

0006  40 

90 

TM  » TMH 

000690 

90 

continue 

OOCuOO 
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Pina 


OPT-l 


FTS  4,b*m4 


o»/i2/^r 


f 0 


ts 


7B 


71 


as 


ta 


bO 


<fi 


101 


103 


IIQ 


113 


120 


CCC3 


100 

UO 


ccs: 


CCS: 


toao 

1 


2 


1 

'» 

10 

11 


12 


1 1 

m 


lODb’} 


H » -DM*<FL0»T(1*1>  - FI 

XJJ  ■ -OH»(FLOATa  - S»J  ♦ II  - FI 

TAIP  • TSEF  - 2.*(rt-riH/10a0. 

SICiA  • SIGUI  » (S1G(J»1I  • SIG(JII*XJJ/10JO. 

P - FPESSUI  * (PKESSlJ^^tl  - FMSSUI)*XJJ/1D00. 

TEHP  • (1TM-0H1»«TEHP*0TE!1PI  ♦ OM'fHII /TH 
aUH  ■ lez,  OlS/TAIRHCl/TEHPI  l/lCt/TEHPI 
IL  • Tm*RLH 

3AS  THE  OROGUE  CHUTE  BEEN  ftCLCASEOT  XF  NOT,  GO  T3  100. 

IF  (ti.lt.tlrcli  so  to  lOO 

HS  ■ MSRGF  - TL  • HO 

OAO  ■ 0i 

SO  TO  110 

HS  ■ HSREF  - TL 

CONTINUE 

TOL  ■ (TM»TEMP|  / (C1»PI 

SUB*  l(6i/PXI*V0LI**  (1./3.I 

OAO  r (PI/4.I*QIAB**2  ».B 

OAS  ■ CAB  * OAO  « SAN 

VE  • •SaRT<ABSIMS/<.SQliaS*SICHA*DASIII 

VEAV  « (VC  * VOI/Ei 

or  * f»oh/veav 

T « T ♦ DT 
VO  ■ VE 

Q ■ .0ail8S»SlGKA*VG**2 
38  • OAB*Q  S to  > OAO*Q  b OH  • OAH»Q 
OS  ■ 08  t 00  * ON 
HM  • HS  ' OS 

HRtTE  ONE  LINE  OF  DATA. 

HRirC  (li3|M(SIGHA|P,TAIR«PI.N,  OT,  Ti  VEt  Qi  TH,  TEHP 
MRXTC  12,41  H,TC,HS,VOLiOXAB,OAS»OAa,OB,DO,DN,  HH 

(F  calculations  ARE  DONE,  GO  BASK  TO  HEAD  HORE  DATA. 

IF  (IFLAG.EOtll  6U  TO  SO 

CONTINUE 


000680 

oouaso 

000700 

aosTio 

000720 
000730 
0007  40 
OC07B1 
000760 
000770 
000780 
CDC7S0 
000800 
000010 
000820 
0038  30 
OQ0040 
000810 
003060 
000870 
OOOOBO 
000800 
oooono 

000910 

ooooto 

000930 
304943 
0009*0 
000960 
303973 
000980 
000990 
001000 
001010 
001020 
0010  33 
001040 
OOIDSO 
3V1063 


001070 

FORNAT  (IHO,*  height  SXGHA(AIN)  PRESSIATHI  TtAlRI  »,  001080 

17X,*L/h»8X*0(TIHEI  TXHE  VEL  (ECU  IL)  *7K  ,1M3 ,7X  ,*‘MASS  (HE  I » , 001090 

2 • TE>1P(HEI»I  u011»3 

FORNAT  I1HQ,«  HEIGHT  TOTAL  LIFT  HT(SVI)  VOL(BALLI»  OOlllO 

1*  0XA(3ALLI  (COSOI  STS  (COSOIBALL  ORACIBALLI  * OOllBO 

ZourAGTCROGI  0RAS(HAXNI  HS  • OS*l  001130 

FORHAT  (1HQ,F7.0,F14.B,F1B.4,  F12.  1,3F1  2.  3,F12.  2,  F 12 . 4,  F12.3  i F U,  1)  3 01140 
rORHAT  ( IHO,  F7 . 0,4K  , 2F 12.  2 ,F  12.3 ,6F12 1 2,  F12.SI  OOllSO 

FORHAT  (IMlf'ALBS  BALLOON  INFLATION  OALOUUATXONS*  , / , 001160 

1 *OPP.OCRAH  P14e*,|X,A10,7FI  001170 

FORHAT  (1H0,11X,*TSHP(HEI  FROM  OEHER  *,3H  ■ ,F7.1,24X,  001100 

1 ’TENPUIRI  • *,F7.1,7,  001190 

2 IHO, lUX, 'TRANSFER  RATE'8X,2Hp  ,FH.  4,  20X,  001200 

5 'INITIAL  PtTIMEI  • ',F9.I»/77I  001210 

format  (1H1,////77,19X,'0RAG  AREA  MAIN  ClUITE  > »,F9.2i  001223 

1 2PX,*URAG  AREA  OROGUE  ■ ',F9.2,/,  301233 


2 1H0,14X,'HT  OPOCUE  AND  HOHRE  ■ ',F7.  3,  001240 
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STOP  301280 

end  001290 
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Nomenclature 


Program  14D 
Cl 

DAB 

DAD 

DAM 

DAS 

DB 

DD 

DH 

DIAB 

DM 

DS 

DT 

DTEMP 


P 

H 

HI 

IFLAG 


■3,0413.  used  in  equation  of  state 
(English  units) 

Drag  area  of  the  balloon,  j 

Drag  area  of  the  drogue, 

Drag  area  of  the  main  chute, 

Drag  area  of  the  whole  system,  (^pSo)syjtem 

Drag  of  the  balloon,  Dg 

Drag  of  the  drogue,  Dq 

Increment  of  height,  AH 

Diameter  of  the  balloon,  d^^ 

Mass  of  helium  transferred  In  time,  At 

Drag  of  the  system,  Dg 

Increment  of  time.  At,  during  which 
system  drops  a height  AH 

Increase  in  temperature,  AT,  of  helium 
due  to  adiabatic  heating  from  the 
balloon's  dropping  a height  AH  (700  ft), 
DTEMP-.8V200  ft  or  4®/1000  ft 

Interpolation  factor  used  to  calculate 
last  line  of  output 

Height,  H 

Initial  height,  H^ 

Flag  used  to  indicate  whether  calcula* 
tions  have  been  completed 
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NOMENCLATURE 


I 

n 

IT0P 

J 

JJ 

P 

PI 

PRESS 

Q 

RLM 

SIG 

SIGMA 

T 

TAIR 

TEMP 

THE 

TL 


Number  of  intervals  of  height  from 
ground , i . e. , I • H/AH 

Index  for  Dfh  loop 

Number  of  height  intervals  of  the 
initial  height  from  the  ground,  i.e., 

IT0P  ■ Hj/AH 

Number  of  lieight  intervals  from 
ground  in  thousands  of  feet. 

Number  of  height  intervals  from  J to  I 

Pressure  (P)  at  height  H 

TT 

Arra/  containing  values  of  pressure  iit 
1000  ft.  intervals,  used  In  the  calculation 
of  P 

2 

Dynamic  pressure,  q ■ >)pv 

Lift  to  mass  ratio  CL/M)  of  helium  In 
the  balloon. 

Array  containing  values  of  o > p/p^ 
at  1000  ft.  Intervals,  used  in  the 
calculation  of  SIGMA 

Ratio  of  0 « p/Pq  height  li 

Time,  t 

Temperature  of  the  air  at  height  H 

Temperature  of  the  helium  in  the 
balloon  at  height  M 

Temperature  of  helium  from  the  dewer 

Total  lift  at  height  H 


37 


TLRHL 


Total  lift  when  drogue  Is  to  be 
released. 


TM 

Total  mass  at  height  H 

TMIIG 

Total  mass  of  helium  to  bo  transferred 
to  balloon 

TMM 

Temporary  storage  for  TM 

TR 

Transfer  rate  of  helium 

TREF 

Temperature  of  air  at 

VE 

Equilibrium  velocity  (v^(H^))  of  system 

VGAV 

Average  velocity  over  Interval 

V0L 

Volume  of  the  balloon 

VO 

Velocity  at  end  of  preceding  Interval 

WD 

Weight  of  the  drogue  and  associated 
hardware 

WS 

Weight  of  the  system  less  the  total  lift 

IVSREF 

Initial  weight  of  the  system 

WW 

■ (WS  < OS)  and  should  always  be  very 
close  to  zero 

XJJ 

Difference  In  height  from  H to  next 
lowest  1000-ft  level 

Adiabatic  Heating  of  the  Gas  in  the  Balloon 


As  the  balloon  descends  into  a thicker  and  thicker 
atmosphere,  the  balloon  contracts  due  to  the  greater  pressure 
causing  the  gas  in  the  balloon  to  undergo  a certain  amount  of 
adiabatic  heating.  This  is  known  to  be  4‘K  per  1000  ft  of 
descent  and  is  essentially  constant  throughout  the  region 
of  interest. 

In  program  P14B  calculations  done  every  200  ft,  and 
through  this  distance  the  gas  (already  in  the  balloon,  not 
in  the  dewer)  will  heat  up  by  .8*K.  The  temperature  and 
transfer  rate  of  the  gas  entering  the  balloon  are  assumed 
to  be  constant,  and  the  gas  in  the  balloon  is  assumed  to  be 
at  all  times  in  equilibrium,  i.e.,  the  temperature  of  the  gas 
is  uniform  throughout  the  balloon.  The  relationship 
connecting  the  temperature,  pressure,  and  volume  of  the 
gas  is  taken  to  be  the  equation  of  state  for  an  ideal  gas 
expressed  as?  _ 

V - 2Ll 

3.0413  P 


where  V is  in  cu  ft,  m in  grams,  T in  ®K,  and  P in 
atmospheres . 

In  each  200-ft  interval  the  system  drops  a height  of 
200  ft  in  a time  6t,  which  is  computed  elsewhere  in  the 
program.  During  this  time  a mass  6m  of  helium  at  temperature 
T,,  is  transferred  at  a uniform  transfer  rate  r,  and  we  have: 


.Wi.i 
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Also,  the  mass  of  helium  already  in  the  balloon  is 
heated  .8“  from  its  original  temperature  to  a temperature 
of  + .8*.  The  mass  m of  helium  entering  the  balloon  is 
at  temperature  and  mixes  to  form  a new  temperature  T 
inside  the  balloon  such  that: 


">o(T  4 .8)  4 AmT„^ 
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II.  VOR  Positioning  Program 


A.  Introduction 

The  position  of  a free  balloon  after  launch  is  found 
by  tracking  the  balloon  from  two  or  more  VOR  stations  which 
measure  the  direction  of  the  balloon  with  respect  to  magnetic 
north.  The  position  of  the  balloon  is  then  found  by  plotting 
these  bearings  on  a map  of  that  region. 

A program  has  been  written  for  the  11P-9810A 
programmable  calculator  which  will  calculate  the  balloon 
position  from  the  station  location.^  and  their  bearings  on 
the  balloon,  Calculations  are  based  on  spherical  geometry, 
which  is  exact  (for  a perfectly  spherical  earth),  and, 
therefore,  errors  due  to  the  plotting  on  the  map  and  to 
distortions  in  the  projection  of  the  map  are  eliminated. 

The  mathematics  are  described  in  Section  B, 

The  program  stores  the  coordinates  of  as  many  as  B 
VOR  stations  at  one  time.  In  this  way  the  path  of  the 
balloon  can  be  followed  over  a longer  distance  without 
having  to  change  the  station  coordinates.  The  program 
works  as  follows: 

After  the  latitudes,  longitudes,  nnd  declinations  of 
the  8 stations  have  been  input,  the  program  will  then 
request  magnetic  bearings  from  any  throe  of  those  stations. 
Tho  user  will  enter  the  number  of  each  station  he  has 
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cKoson  iiiul  itH  corrosponding  boaring.  Taking  tlio  .stations 
two  at  a time,  the  program  will  then  calculate  intorsoctlon 
points  among  the  throe  station.^  scloctod. 

After  the  coordinates  of  these  intersection  points 
l\ave  been  calculated  and  printed,  the  program  will  request 
a now  set  of  bearings.  These  new  bearings  need  not 
necessarily  be  from  the  same  stations  previously  chosen; 
any  3 of  the  6 may  be  selected. 

A complete  description  on  how  to  use  the  program  is 
given  in  Section  C. 


B.  Mathematics 
1.  Introduction 

The  projection  on  the  earth  of  the  line  of  sight  from 
a VOR  station  to  the  balloon  follows  a great-circle  path. 
Therefore,  the  balloon  position  can  be  found  by  calculating 
the  location  of  the  point  of  intersection  of  the  two  great 
circles  which  correspond  to  the  lines  of  sight  from  fixes 
from  the  two  given  VOR  stations. 

Inputs  to  the  calculation  are  the  latitudes,  longitudes, 
and  declinations  of  the  two  stations  and  the  directions  of 
the  balloon  from  the  stations  with  respect  to  magnetic 
north.  The  declination  is  used  to  change  the  magnetic 
direction  to  a compass  heading  with  respect  to  true  north. 

For  each  station,  these  quantities  are  used  to  form 
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ono  vector  which  points  from  the  station  to  the  balloon, 
and  another  vector  which  points  from  the  center  of  the 
earth  to  the  station.  Taking  the  cross  product  of  these 
two  vectors  produces  a vector  which  is  perpendicular  to 
the  plane  of  that  great  circle  which  is  the  projection 
of  the  line  of  sight  from  that  station  to  the  balloon.  In 
other  words,  one  set  of  station  coordinates  and  a bearing 
from  that  station  produces  a vector  which  is  normal  to 
the  plane  of  that  great  circle  which  corresponds  to  the  line 
of  sight  from  that  station  to  the  balloon. 

Taking  the  cross  product  of  two  of  these  normal 
vectors  produces  a vector  which  points  from  the  center  of 
the  earth  towards  the  point  of  intersection  of  the  corres- 
ponding two  great  circles,  l.e.,  it  points  along  the  line 
of  Intersection  of  the  two  great-circle  planes.  From  this 
vector  one  can  compute  the  latitude  and  longitude  of  the 
balloon . 

2,  Input  and  Output 

The  input  to  the  calculation  is  the  latitudes,  longi- 
tudes, declinations,  and  magnetic  bearings  of  the  balloon 
of  two  stations. 

The  latitudes  are  north  latitudes. 

The  longitudes  are  input  to  and  output  from  the  program 
as  west  longitudes.  However,  for  mathematical  convenience, 
within  the  program  and  in  the  calculations  which  follow, 
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east  longitude  is  used  In  order  to  be  consistent  with  a 
right-handed  Cartesian  coordinate  system. 

In  the  program  the  magnetic  bearings  are  converted  to 
bearings  with  respect  to  true  north  by  adding  to  each  its 
corresponding  declination  (expressed  as  positive  for  an 
easterly  declination). 

The  north  latitude  and  west  longitude  of  the  balloon 
position  are  printed  by  the  program. 


3.  Nomenclature 


0 

A,  C 
B 

®A'®B’®C 

^A’^B'^C 
OA,  OB,  0^ 


clJ 


AE.  CE 


is  the  center  of  the  earth 

are  two  VOR  stations 

is  the  position  of  the  balloon 

are  the  north  latitudes  of  positions  A,  B, 
and  C,  respectively 

are  the  east  longitudes  of  positions  A,  B, 
and  C,  respectively 

are  unit  vectors  from  0 to  A,  B,  and  C, 
respectively.  These  vectors  are  local 
verticals 

are  unit  vectors  at  A and  B,  respectively, 
which  point  due  north,  l.e.,  these  vectors 
are  local  horizontals  pointing  north 

are  unit  vectors  at  A and  B,  respectively, 
which  point  due  east,  l.e.,  local 
horizontals  pointing  east 
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n^iHg  are  bearings  on  the  balloon  with  respect  tc  true 

north  from  stations  A and  C,  respectively 

At,  Ct  are  unit  vectors  which  point  from  A and  C, 

respectively,  along  those  great-circle  paths 
which  are  in  the  direction  of  the  balloon  B, 
l.e.,  local  horizontals  pointing  in  the  direction 
of  the  balloon 

are  vectors  which  are  normal  to  the  great-circle 
planes  specified  by  the  vectors  ^ and  CB, 
respectively 

4.  Equations 

a.  The  first  step  is  to  form  the  vectors  oit, 

A^,  and  for  station  A.  The  vectors  jyt!,  and  will 
then  be  constructed  in  a similar  fashion.  The  vectors 
are  expressed  in  terms  of  their  x,  y,  and  z components 
respectively.  The  origin  of  this  Cartesian  coordinate 
system  lies  at  the  center  of  the  earth;  the  z-axls  points 
towards  the  north  pole;  the  x and  y axes  lie  jn  the  plane 
of  the  equator  and  point  toward  the  0^  and  90‘  meridians, 
respoctlvely . 

The  X,.  y,  and  t components  of  these  vectors 
can  be  expressed  in  terms  of  the  latitude  6^  and  longitude 
of  station  A. 

Looking  at  Figures  1 and  2,  one  can  see  that  the 
X,  y,  and  z components,  respectively,  of  these  vectors  are: 

0^  ■ (co80^  cosX^,  cos®^  sinX^,  slne^  ) 
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a15  » (-sin6y^  coaX^,  -alne^  ainX^,  cos6^) 

■ (-sinX^,  cosX^,  0) 

The  vectors  0?,  Cl5,  and  cS  are  constructed 
simi 1 iarly , and  their  components  can  be  found  by  substituting 

^A  vectors  given  above. 

b.  The  next  step  is  to  form  the  vector  ^ which 
points  from  A along  the  great-circle  path  which  Is  the 
projection  of  the  line  of  sight  from  A to  B.  In  Figure  3 
is  a diagram  of  the  local  horizontal  plane  at  A.  From  this 
we  see  that: 

a5  ■ ANf  cosn^^  + A?  sinn^ 

Expanding  this  to  find  the  x,  y,  and  z components  of  we 
hove ! 

■ C-sinOy^  cosX^  coan^  - sinX^  sinn^, 

-sinQ^  sinX^  cosn^  cosX^  Binn^, 
cose^  costly^  ) 

The  vector  ^ is  constructed  similiarly,  and  its 
components  can  be  found  by  substituting  0^,  for  0^,  X^,  for  X^ 
and  for  in  the  above  equations. 

c.  The  next  step  is  to  form  the  vector  n*^  which  is 
normal  to  that  groat-circle  plane  which  contains  the  lino  of 
sight  from  station  A to  balloon  B.  This  vector  can  be  formed 
by  crossing  uny  two  vectors  which  lio  in  the  plane  (and  are 
not  purallolj.  The  vectors  5\  and  Alt  lie  in  this  plane,  and 
wo  have: 


I 


I 


1 


•1 


•1 


I 


i 
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"a  ■ 

X JTt 

(n“) 

■ 

X 

(CT),  - 

(At), 

■ 

slnX^ 

COSHy^  - 

•*”®A 

cosX 

(n-;) 

■ 

y 

(rt)x  ■ 

(5J), 

(A») 

■ 

-COSX^ 

1 “*"a 

- sinBy 

^ sin 

Cn-J) 

■ 

t 

Cffi«)x 

(Tt)y  • 

(«)|, 

m) 

m 

co.Sa 

“""A 

The  vector  is  conetrueted  slmlllerl/.  end  Ite 
component!  can  be  found  by  substituting  6^  for  6^.  for  X^ 
end  for  n^  in  the  above  equations. 

d.  Taking  the  cross  product  ^ of  the  vectors  n^  and 

n^  produces  a vector  which  points  along  the  line  of  Intersection 
of  the  two  great-circle  planes.  This  line  of  Intersection 
passes  through  the  center  of  the  earth  and  the  balloon.  The 
vector  formed  will  point  either  from  the  center  of  the  earth 
to  the  balloon  towards  a north  latitude  or  In  the  opposite 
direction  towards  a south  latitude.  The  components  of  are; 

% • <Vy 

’y  ■ '"b>. 

'’l  ' '“*B)y  • '"aJ  '^B>x 

In  the  ease  that  7 points  towards  a south  latitude, 
we  can  turn  the  vector  around  by  making  the  transformation 
V ■>  -V. 

e.  Even  though  v now  points  In  the  direction  of  OB, 


V must  bo  divided  b/  Its  length  and  we  have: 

ol'  ■ V / 1^1 

whore  |v|  ■ + v‘ 

X y z 

6)  Looking  back  at  Figures  1 and  2,  the  coordinates 
of  tho  balloon  3^  and  can  bo  found  from  0^  or  7 by; 


B 


'B 


- sin‘^(ff«)^]  . sln"^  j^) 

1.  .1 


Ag  - cos 


[ 


cos6, 


•]  ■ cos  [ 


|7| COS9, 


C.  The  Program 

1 . Introduction 

The  following  is  a step-by-step  description  of  how 
to  use  the  VOR  positioning  program.  A sample  output,  a 
listing,  a flow  chart,  and  an  explanation  of  the  organization 
of  the  memory  are  also  provided  in  Figures  4-7,  respectively 
Supplementary  information  about  the  program  is  given  which 
enables  the  user,  when  data  has  been  entered  incorrectly,  to 
correct  the  mistake  without  having  to  go  back  to  tho  beginning 
of  the  program. 

2.  How  to  Use  the  Program 

Press  END,  CONTINUE. 

The  program  then  asks  for  the  number  of  stations 
to  be  put  into  tho  calculator.  Tho  program  can  handle  up  to 
8 stations.  Tho  user  puts  this  number  into  the  x register 
and  then  presses  CONTINUE. 

Tho  program  then  prints  a station  number,  starting 
at  station  1 and  going  up  to  the  number  of  stations  requested 
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In  th6  previous  step,  and  requests  the  latitude,  longitude, 
and  declination,  respectively,  of  that  station.  For  each 
of  these  variables,  enter  the  number  of  degrees  in  z, 
minutes  in  y,  and  seconds  in  x,  and  then  press  CONTINUE. 

Latitude  is  input  as  north  latitude,  longitude 
as  west  longitude,  and  declination  is  input  as  a positive 
quantity  for  an  easterly  declination  and  as  a negative 
quantity  for  a westerly  declination. 

After  reading  these  quantities  for  a particular 
station,  the  program  goes  on  to  the  next  station,  or,  if  the 
last  station  has  just  been  Input,  it  goes  on  to  the  next  step. 

The  program  then  asks  for  the  time  of  day.  This 
is  done  so  that  there  will  be  a record  of  the  time  that  the 
VOR  readings  were  taken.  The  user  puts  the  time  Into  the 
X register  in  the  following  way:  14:45  is  put  in  as  14. 4S, 

etc.  The  user  then  presses  CONTINUE. 

The  program  then  asks  for  bearings  from  any  three 
of  the  stations,  called  1ST,  2ND,  and  3RD  stations,  respectively 
For  each  of  the  three  stations  to  be  selected,  enter  the 
number  of  the  station  required  in  the  y register  and  its 
corresponding  magnetic  bearing,  expressed  in  decimal  degrees, 
in  the  X register. 

The  program  then  prints  the  coordinates  of  the 
balloon  location  calculated  using  the  1st  and  2nd  stations, 
then  the  1st  and  3rd  stations,  and  finally  the  2nd  and  3rd 
stations.  For  each  station  pair  the  latitude  and  longitude 
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are  printed  In  degrees,  minutes,  and  seconds. 

The  program  then  asks  for  a new  time  and  then 
new  station  numbers,  and  bearings.  Again,  any  throe  stations 
may  be  selected,  and  these  need  not  be  the  same  ones  which 
were  solectod  previously.  After  each  set  of  bearings  has 
been  Input,  and  the  balloon  positions  calculated,  the 
program  goes  back  to  ask  for  a new  tine  and  a new  set  of 
station  numbers  and  bearings. 

A sample  output  is  shown  In  Figure  4.  All 
variables  which  the  user  must  supply  have  an  asterisk 
placed  beside  them.  Note  that  all  variables  input  by  the 
user  are  printed  on  the  paper  tape  Immediately  after 
being  entered. 

3.  Supplementary  Information 

In  order  to  help  the  user  understand  (and  j 

possibly  modify)  the  program,  a listing  of  the  program  is 

given  In  Figure  5 and  a flow  chart  in  Figure  6.  ' 

In  the  case  that,  after  the  station  data  have 
boon  entered,  it  la  discovered  that  a station  coordinate 
was  ontorod  incorrectly,  or  it  is  necessary  that,  another 
station  be  substituted  or  addod,  it  is  usually  more 
convenient  to  bo  able  to  enter  the  information  directly 
Into  the  memory  rather  than  having  the  program  do  it. 

(For  the  program  to  do  this  requires  that  all  of  the 
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Information  b«  roenterdd . ) This  requires  a better 
knowledge  of  the  memory  and  program  organization.  The 
organization  of  the  memory  is  shown  in  Figure  7. 

To  put  a number  into  a particular  memory 
location,  put  the  number  into  the  x register,  press  x->’CXTO) 
followed  by  the  three-digit  memory  location.  Since  the 
coordinates  entered  in  this  program  are  always  expressed 
in  terms  of  degrees,  minutes,  and  seconds,  the  user  must 
first  change  this  to  dsclmal  degrees  before  entering  it 
into  memory.  In  this  case,  put  the  numbers  of  degrees, 
minutes,  and  seconds  into  the  z,  y,  and  x registers, 
respectively,  press  K,  6 to  convert  to  decimal  degrees  in 
the  x register,  and  then  press  x-^,  followed  by  the  three- 
digit  memory  address. 

For  station  number  i,  i ■ 1 to  6,  the  latitude  e, 
longitude  X,  and  declination  n are  put  into  memory  locations 
Si  + 5,  5i>6,  and  5i'«’7,  respectively.  The  number  of  stations 
which  are  being  stored  in  the  memory  is  kept  in  memory 
address  000. 

If  a bearing  is  entered  Incorrectly,  press  GTO, 

0.  2,  5,  0,  and  CONTINUE.  The  program  will  then  ask  for  a 
time  followed  by  three  new  station  numbers  and  their 
corresponding  bearings.  The  old  station  numbers  and  bearing 
are  rt^placod. 
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D,  Comparison  of  VOR  Data  with  Radar  Data 

Using  data  from  a particular  balloon  flight  for  which 
both  VOR  and  radar  data  were  available,  a comparison  of 
results  using  these  two  sots  of  data  was  made.  Radar  data 
are  considered  to  produce  substantially  more  accurate 
results  than  VOR  data  and,  therefore,  are  suitable  for 
determining  the  accuracy  of  the  VOR  data. 

Errors  in  VOR  positioning  can  result  from  a number  of 
sources:  1)  the  error  (An)  in  the  magnetic  bearing  (n) 

of  the  target  as  viewed  by  the  VOR  station  in  question, 

2)  interference  from  other  radio  signals  or  mountains,  and 

3)  the  line  of  sight  being  too  close  to  the  horizon.  The 
error  in  the  bearing  is  the  only  type  of  error  which  will 
be  discussed  here.  Errors  due  to  the  other  two  sources 
are  highly  unpredictable  and  would  probably  have  to  be 
measured  empirically. 

The  error  in  positioning  (E)  due  to  one  VOR  station 
alone  at  a distance  d from  the  target,  assuming  the  bearing 
of  the  other  station  used  for  the  fix  is  exact,  is: 

E ■ d tan(An),  or  for  small  An,  E <=  d(An).  The  maximum 
error  in  a VOR  bearing  is  known  to  be  +_  112^.  Errors 
(iT^  and  1^)  , due  to  the  errors  in  the  bearings  of  both  of 
the  VOR  stations  used  to  produce  a given  fix,  add  vectorially 
to  produce  the  total  error  F = E*j  + E*, . This  is  shown 
in  Figure  8. 


53 


- ' ujeMiiltt rtlMiiinii 


When  the  angle  a,  the  angle  between  the  lines  of  sight 
from  two  given  VOR  stations,  approaches  180*,  i.e.,  when  the 
balloon  is  very  close  to  the  line  Joining  the  two  VOR 
stations,  the  total  error  can  become  very  large 

even  though  the  bearings  themselves  are  accurate,  i.e., 

An^  and  Ari2  are  small.  Conversely,  an  intersection  angle 
of  90*  between  the  lines  of  sight  is  optimum.  This  is  shown 
in  Figure  9.  The  shaded  areas  show  where  the  balloon  might 
be  for  each  of  the  two  fixes  A and  B,  assuming  an  error  of 
^ 1/2*  in  each  of  the  bearing-angle  measurements.  The 
centers  of  the  shaded  areas  are  the  positions  where  the  VOR 
predicts  that  the  balloon  is  located.  Note  that  for  position 
B the  fact  that  o Is  close  to  180*  can  lead  to  a large  error 
in  the  prediction  of  the  balloon  position. 

Both  radar  and  VOR  data  were  available  on  balloon  flight 
H75-4,  which  took  place  on  February  7,  1975.  VOR  data  were 
available  from  the  following  stations:  1)  Truth  or  Consequences, 

2]  Doming,  3)  Newman,  Texas,  4)  Holloman,  5)  Pinon, 

6)  Hudspeth,  Texas,  7)  Roswell,  and  8)  Pecos,  Texas.  Those 
stations  and  the  position  of  the  radar  are  shown  on  a map  of 
the  area  in  Figure  10.  Balloon  positions  and  their  corresponding 
times  are  shown  throughout  the  flight. 

During  this  flight  balloon  positions  were  computed  using 
Truth  or  Consequences,  Holloman,  Pinon,  and  Roswell.  The 
other  stations  were  too  far  away  to  provide  accurate  results. 
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Since  B > d tan  An,  the  error  Is  proportional  to  the  distance 
of  the  station  from  the  target.  Furthermore,  the  farther 
away  the  station  is  from  the  balloon,  the  closer  to  the 
horizon  the  balloon  appears. 

Larger  discrepancies  ( 2 to  4 nm]  occurred  at  the 
beginning  of  the  flight  (before  9:14)  when  the  balloon  had 
not  yet  reached  flight  altitude  and  was  therefore  too  close 
to  the  horizon  to  get  a good  fix.  The  radar  was  much 
closer  to  the  balloon  during  this  time,  and  its  results  are 
probably  accurate. 

Later  in  the  flight  (10:35  and  after),  discrepancies  of 

2 to  5 nm  were  seen.  At  this  point  the  balloon  was  about 
h.'\lfway  between  Roswell  and  Plnon,  so  that  o was  close  to 
ISO*.  Therefore,  these  two  stations  could  not  be  used  together 
to  get  a good  fix.  Holloman  and  Truth  or  Consequences  were 

at  this  time  quite  far  away  (60  and  120  nm,  respectively) 
and  were  also  looking  over  a mountain  range.  Hence,  using 
either  one  of  these  with  either  Roswell  or  Pinon  also  tended 
to  produce  larger  errors  than  usual  because  of  both  the 
distance,  ircreasing  the  error  and  making  the  balloon 
appear  close  to  the  horizon,  and  the  Interference  produced 
by  the  mountains.  This  seems  to  have  boon  in  fact  the  case, 
because  the  fixes  using  Holloman  were  consistently  about 

3 nm  to  the  south,  whereas  the  fixes  using  Troth  or 
Consequences  were  consistently  about  3 nm  to  the  north.  Note 
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that  the  radar  data  could  alao  have  been  Influenced  by  the 
mountains  and  the  balloon's  being  too  close  to  the  horizon. 
However,  the  fact  that  Holloman  predicted  positions 
were  consistently  6 nm  south  of  those  predicted  by  Truth 
or  Consequences  shows  that  the  VOR  data  themselves  were 
inconsistent.  During  this  part  of  the  flight  Carlsbad 
could  have  provided  a good  fix,  but  data  from  it  were  not 
available . 
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in.  Totherod- Bal  loon  Experiment 
Introduction 

A series  of  testa  was  made  in  Februarx-March  1975  to 
aid  in  the  implementation  of  the'Tactlcal  Loran  Tritether 
Instrumentation  Hook-up',  These  tests  consisted  of  flying 
a 30000  cu  ft  balloon  on  tethers  of  various  lengths  ClOO, 

500,  1000,  and  2000  ft)  with  various  loads  (from  no  load 
up  to  260  lbs)  in  various  winds.  The  purpose  of  these  tests 
was  to  determine  the  extent  of  balloon  wander  under  different 
winci  conditions,  The  balloon  position  was  measured  by 
moans  of  theodolites,  and  wind  azimuth  and  velocity  were 
measured  by  Instruments  attached  just  under  the  balloon. 

Table  1 gives  a summary  of  the  different  runs  made  In  this 
series  of  tests. 

Analysis 

At  the  beginning  the  approach  was  to  plot  Individual 
data  points,  with  each  graph  containing  data  from  only 
one  run.  Balloon  positions  with  respect  to  ground  zero 
were  plotted  with  the  y-axls  representing  north  and  the 
x-uxls  oast.  Values  of  wind  direction  and  azimuth 
corrospondlng  to  each  balloon  position  were  also  plotted  on 
thosn  graphs,  with  the  distance  from  tho  origin  corresponding 
to  tho  wind  velocity  and  tho  angle  with  respect  to  the 
y-nxln  correspondinR  to  tho  wind  aaimutl, . The  two  sots 
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of  points  were  numborcJ  and  connected  by  straight  linos, 
providing  a kind  of  time  history  of  the  balloon  position, 
wind  direction,  and  wind  speed. 

Using  the  corresponding  values  of  wind  azimuth,  the 
balloon  positions  were  then  resolved  into  components  parallel 
and  perpendicular  to  the  wind.  These  values  were  also 
plotted  as  a function  of  time. 

Since  not  much  could  be  inferred  from  those  types  of 
plots,  since  the  wind  data  might  not  have  been  good  enough  to 
justify  a point-to-point  breakdown,  and  since  it  takes  a 
certain  amount  of  time  for  a given  change  in  wind  direction 
or  speed  to  effect  a change  in  balloon  position,  it  was 
decided  to  group  together  those  times  when  both  wind  speed 
and  direction  were  'approximately  constant.'  Balloon  positions 
durinu  a given  time  period  were  resolved  Into  components 
parallel  and  perpendicular  to  the  wind  baaed  on  the  average 
azimuth  of  tho  wind  over  the  whole  period.  Such  time  period 
contains  at  least  five  points,  and,  for  some  runs  in  which 
the  wind  was  nearly  constant,  it  was  not  necessary  to  divide 
tho  data  ut  all.  Balloon  positions  parallel  und  perpendicular 
to  the  wind  wore  graphed  versus  time,  and  each  run  was 
graphed  separately. 

Sinco  this  type  of  tiiiio  history  of  balloon  wander  did 
not  seem  to  provide  any  notubU'  insights,  data  from  each 
of  those  time  periods  wore  grouped  logothor,  and  statistics 
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of  mean,  standard  deviation,  maximum  excursion  from  ground 
zero,  and  spread  of  excursions  about  ground  zero  were 
calculated.  These  are  statistics  of  the  balloon  displacement 
from  ground  zero  both  parallel  and  perpendicular  to  the  wind. 
Hence,  for  each  time  period  of  a given  run  one  value  for 
each  of  the  above  statistics  of  balloon  displacement  was 
calculated,  both  parallel  and  perpendicular  to  the  wind. 

In  the  next  set  of  graphs,  each  of  the  statistics  of 
displacement  listed  above  was  plotted  versus  the  average  wind 
velocity  of  the  corresponding  time  interval.  Bach  graph 
contained  the  data  from  all  runs  of  a given  tether  length, 
for  a given  statistical  parameter  either  parallel  or 
perpendicular  to  the  wind.  If  a particular  run  was  divided 
into  more  than  one  time  interval,  more  than  one  set  of 
statistics  were  plotted  for  that  run. 

Note  that  in  the  case  of  maximum  and  minimum  excursions, 
the  length  of  the  time  interval  of  the  run  is  important; 
however,  no  attempt  was  made  to  take  this  into  consideration. 
On  the  other  hand,  in  the  case  of  the  mean  and  standard 
deviation,  the  length  of  the  time  interval  is  not  important 
as  long  as  it  Is  'long  enough*  to  provide  a representative 
sample . 

Prom  this  set  of  graphs  it  was  seen  that  in  the  runs  of 
the  100'  tethor  the  balloon  wander  is  a much  higher  fraction 
of  the  cnblo  length  than  with  the  other  tether  lengths. 
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Noto  thntj  in  the  application  for  which  this  will  be  used,  it 
is  the  conith  angle  of  the  balloon  tether  which  Is  important. 

Fr  om  this  sot  of  graphs  it  was  clocidotl  that  the  100'  tether 
was  not  suitable, 

One  can  also  see  that  tho  1000*  and  2000*  tethers  are 
noticeably  bettor  than  the  500'  tether.  However,  it  was 
felt  that  performance  of  these  lengths  would  not  justify 
their  iso  over  the  500'  tether,  Use  of  u longer  cable 
increases  tho  weight,  thereby  decreasing  the  net  lift  of  tho 
balloon,  und  makes  the  rigging  of  tho  balloon  more  ill  ffi  cult 

Since  it  is  tho  zenith  angle  of  tho  balloon  tothor  which 
Is  important  in  this  appllcaiton,  tho  statistics  of  balloon 
displacement  listed  above  wero  converted  to  percentage 
displacement  (■dlsplacoinont/cuble  length  x 100).  Each 
statis  ic  was  plotted  soparaiely  versus  average  wind  velocity. 
For  a given  statistic,  tho  data  for  the  500',  1000',  and 
20110'  tothors  parallel  or  perpendicular  to  the  wind  were 
plotted  on  tlio  same  graph. 

For  moan  d i .s p 1 nc ement  parallel  to  tho  wind,  the  poreenttj’, o 
d i li  |)  I at  .MiiL’n  t inoroasos  with  wind  velocity  In  a way  which  la 
almost  linear.  Tliere  is  no  noticoablc  dcpoiulenco  on  tether 
loniith  apart  From  tho  fact  that  the  500'  tether  doe.s  not 
perform  c|Uite  us  wo’l  as  tho  two  longer  tethers,  Tliorc  is 
a fair  amount  of  scatter  in  tho  graph,  but  tlio  slope  of  that 
lino  which  produces  tho  bo.st  fit  is  1.7%  of  cable  length  per 


mph  of  average  wind  velocity.  Therefore,  the  degree  to 
which  the  balloon  is  blown  back  into  the  wind  is  proportional 
to  the  average  wind  velocity. 

Mean  displacement  perpendicular  to  the  wind  showed  no 
systematic  variation  with  wind  velcolty.  This  is  reasonable 
since  the  balloon  should  not  have  a preferred  average 
position  perpendicular  to  the  wind. 

Standard  devlalton  parallel  and  perpendicular  to  the 
wind  is  very  variable  and,  with  the  exception  of  windS' 
below  S mph,  shows  no  variation  with  wind  velocity.  Further- 
more, the  SOO*  tether  data  have  noticeably  higher  standard 
deviations  than  do  the  1000'-  or  2000'-  tether  data.  Note 
that  even  in  high  winds,  greater  than  20  mph,  the  standard 
deviation  did  not  increase. 

The  maximum  excursion  from  ground  zero  parallel  to  the 
wind  showed  an  almost  linear  Increase  with  wind  velocity, 
but  wi '.h  considerable  spread.  Again  the  500'  tether  shows, 
in  gfinoral,  greater  displacements  than  do  the  two  longer 
tethor  lengths. 

The  spread  of  excursions  parallel  to  the  wind  shows  no 
relation  to  wind  velocity  above  5 mph.  But,  here  again,  the 
500'  tdther  shows  greater  displacements. 

The  variation  of  these  two  statistics  perpendicular  to 
tho  wind  also  showed  no  variation  with  wind  velcolty  above 
5 mph,  but  tho  500'  tothor  showed  significantly  greater 
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displftcoments  than  did  the  two  longer  tether  lengths. 

We  have  seen  that  parallel  to  the  wind  the  mean  and 
nm.xlmum  excursion  varied  with  wind  velocity  in  an  almost 
linear  fashion.  Dividing  each  of  those  statistic*  by 
wind  velocity,  (producing  a mean  or  maximum  displacement 
/mph  average  wind  velocity),  these  data  were  plotted  versus 
load  to  try  to  find  what  effect,  if  any,  the  load  might 
have  on  balloon  wander.  Qualitatively,  we  know  that  an 
increase  in  the  load  will  Increase  the  angle  of  attack  of 
the  balloon.  Up  to  a certain  point  this  is  desirable  'ilnce 
it  produces  an  Increase  in  aerodynamic  life,  but  after  this 
point  the  aerodynamic  lift  derreases  with  angle  of  attack. 

For  the  "mean  displacement  parallel  to  the  wind/mph 
average  wind  velocity",  there  seemed  to  be  some  correlation 
with  load,  i.e.,  angle  of  attack.  The  spread  is  very  large, 
but  one  can  conclude  that  the  slope  of  the  percentage  mean 
displacoment  parallel  to  the  wind  vei'sus  wind  velocity  is 
higher  when  the  load  is  greater. 

Tho  slopes  computed  using  the  percent  ugo-maxlmuir- 
oxcursion  statistics  show,  on  the  other  hond,  no  corrolutiua 
with  loud. 

Summary 

1)  Data  were  first  resolved  into  components  parallel 
and  perpendicular  to  the  wind. 


2)  Data  were  then  grouped  together  into  time  periods 
when  both  tho  speed  and  azimuth  of  the  wind  were  constant. 

3)  Data  for  each  time  period  during  which  the  wind 
was  constant  were  then  grouped  together  to  form  statistics 
of  the  displacement  of  the  balloon  with  respect  to  ground 
zero  both  parallel  and  perpendicular  to  the  wind. 

4)  The  percentage  displacement  of  the  100*  tether  was 
too  high  to  be  suitable  for  this  application.  The  percentage 
displacements  for  the  1000*  and  2000*  tethers  were  noticeably 
smaller  than  those  of  the  500',  but  their  performance  was 

not  that  much  bettor  to  justify  their  use  over  the  shorter 
500  * t^'ther. 

1)  Percentage  moan  displacement  parallel  to  the 
wind  varied  linearly  with  wind  velocity,  with  a slope  of 
1.7^/mph.  the  balloon  being  blown  back  into  the  wind  more 
and  more  as  the  wind  increased.  Perpendicular  to  the  wind 
tho  mean  displacement  showed  no  variation  with  wind  velocit/, 
l.e.,  the  balloon  had  no  preferred  position  perpendicular 
to  thr  wind, 

6)  Maximum  excursion  parallel  to  tho  wind  varies 
linouriy  with  wind  velocity.  This  is  due  to  a great  extent, 
however,  to  the  variation  of  tho  moan.  Perpendicular  to  the 
wind  thoro  Is  no  syitomatic  variation. 

7)  Standard  deviation  purollel  and  perpendicular  to 
tl\o  wind  showed  no  variation  with  the  wind  above  5 mph . 


SO 


8)  In  the  case  of  the  standard  deviation  parallel  and 
perpendicular  to  the  wind  and  of  the  maximum  excursion  parallel 
to  the  wind,  the  500'  tether  showed  greater  percentage  dis- 
placements, in  general,  than  did  the  two  longer  tethers, 

9)  Load,  in  general,  seems  to  produce  greater  displacements 
parallel  to  the  wind. 
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Table  1 


2000  ' 
tether 


1000' 

tether 


500 ' 
tether 


100  ' 
tether 


NO 

LuAD 

2/26  #4 
2/27  #4 
3/4  M* 

2/26  #3 
2/27  #3 
3/4  #3 

2/26  #2 
2/27  #2 
3/4  #2 

2/26  #1 
2/27  #1 
3/4  n 


LOAD 

BAR  ONLV 

3/7  #2* 
3/7  #3 

3/7  #4* 


2/26  #5 
3/3  #5 

3/3  #4 


LOAD  BAR 
+14C  LS5 
BALLAST 

2/27  #8 

2/28  #4 

3/6  #4 

2/27  #7 
2/28  #3 
3/6  #3 

2/27  #6 
2/28  #2 
3/6  #2 

2/27  #5 
2/28  #1 
3/4  #5 

3/6  #1 


LOAD  BAR 
+265  LBS 
BALLAST 

3/7  #1 


3/3  #3 


3/3  #2* 
3/7  #6 

3/3  #1* 
3/7  #5 


* NO  WIND  DATA 
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IV  Salloon-Pllght  Simulation 

The  following  balloon-flight  simulations  were  done. 

A,  Eglin  APB  simulation 

Using  10  years  of  wind  data  from  Bglin  APB,  balloon 
positions  and  statistics  of  balloon  positions  were  calculated 
assuming  an  ascent  to  10  mb  at  a constant  ascent  rate  of  800 
ft/min,  letting  the  balloon  float  for  0,  1,  4,  6,  8,  and  10 
hours  at  10  mb.  This  was  done  using  wind  data  for  days  1- 
10,  11-20,  and  21-31  in  each  of  the  months  from  April  to 
September . 

The  wind  data  were  received  on  a 9-track,  phase-encoded 
tape  which  had  to  be  converted  to  a 7-track  tape  in  order  to 
make  it  compatible  with  the  CDC  6600  at  APGL,  Before  writing 
the  7-track,  800-bpl  tape  the  data  were  sorted  by  month  and 
day . 

The  balloon  positions  and  associated  statistics  were 
done  by  a program  called  CHICSIM,  which  was  modified  for  this 
particular  job. 

B.  }lolloman  and  White  Sands  simulation 

Using  10  years  of  wind  data  from  Holloman  and  White  Sands, 
balloon  positions  and  statistics  of  balloon  positions  were 
calculated  assuming  ascents  to  80,  90,  and  100  kft,  at  a 
constant  ascent  rate  of  800  ft/min,  lotting  the  balloon  float 
for  0,  2,  A,  6,  8,  10,  and  12  hours  at  each  of  those  altitudes. 
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Ench  of  the  nbovo  cases  was  run  for  days  1-15  and  J6-31  for 
all  months  of  the  calendar  year. 

These  wind  data  were  already  on  a CDC  7-track  tape, 

The  data  were  processed  by  a program  called  SIMBALL,  which 
was  modified  slightly  for  this  particular  job. 


